Eccentricity fluctuations and its possible effect on elliptic flow measurements 
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The elliptic flow measured at RHIC has been interpreted as a signature for strong partonic in- 
teractions early in the collision and as an indication of a well developed quark-gluon plasma phase. 
The measured values of elliptic flow, using methods based on multi-particle correlations, are af- 
fected by fluctuations in the magnitude of the elliptic flow. In this Letter, using a Monte Carlo 
Glauber calculation, we estimate what the possible effect of spatial eccentricity fluctuations is on 
the determination of elliptic flow. 

PACS numbers: 25.75.-q, 25.75.Ld, 25.75.Dw, 25.75.Gz, 24.10.Lx 



In non-central heavy- ion collisions, the initial spatial 
anisotropy due to the geometry of the overlap region and 
the pressure developed early in the collision generate an 
observable azimuthal momentum-space anisotropy. The 
particle yields produced in heavy-ion collisions can be 
characterized by 
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where pt is the transverse momentum of the particle, <f> is 
its azimuthal angle, y is the rapidity and the reaction 
plane angle, see fig [J The second coefficient, v 2> of this 
Fourier series is called elliptic flow. 
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FIG. 1: Schematic view of a nucleus- nucleus collision in the 
transverse plane. 



Elliptic flow as a signature of hydrodynamic behavior 
of nuclear matter produced in high energy nuclear colli- 
sions has been proposed by Ollitrault in 1992 After 
that it has attracted increasing attention from both ex- 
perimentalists and theorists [5| and has been measured 
at AGS [E0, SPS HE3 and RHIC [HEDS 

energies. 

It is thought that elliptic flow reflects the amount of in- 
teractions between the constituents at an early time in 
the evolution of the produced system . Therefore it is 



sensitive to the equation of state of the produced system 
when this system might be in the quark-gluon plasma 
phase. 

Since the reaction plane is not known experimentally, 
the elliptic flow is calculated using azimuthal angular cor- 
relations between the observed particles . In the case 
of two particle correlations the measurement is propor- 
tional to v 2 . The reported elliptic flow values are there- 
fore obtained as \J (v 2 ) after averaging over events. 

Because elliptic flow is a collective effect, it is a corre- 
lation of all the particles with the reaction plane. This 
can be exploited experimentally by using multiple parti- 
cle correlations to calculate v 2 . To calculate these cor- 
relations a convenient mathematical approach is to use 
cumulants. This method, proposed in [l£j, has the addi- 
tional advantage that it allows to subtract the so called 
non-flow effects from v 2 . Non-flow effects are correla- 
tions between the particles not related to the reaction 
plane. Such effects include, but are not limited to, reso- 
nance decays, (mini)jet fragmentation and Bose-Einstein 
correlations. The cumulant method uses multi-particle 
correlations which introduce higher powers of v 2 . The 
corresponding equations for calculating v 2 in the cumu- 
lant method for two, four and six particle azimuthal cor- 
relations are given by: 
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= -;((v!)-9(vl)(v 2 2 ) + 12(v 2 2 f), (2) 



where v 2 {2} is calculated using two-particle azimuthal 
correlations, t>2{4} using a mix of two and four-particle 
azimuthal correlations and i>2{6} uses two, four and six- 
particle azimuthal correlations. 

However, due to event-by-event fluctuations in the el- 
liptic flow for instance, the event averaged (v 2 ) ^ {v 2 ) n 
for 71 > 2 17]. Therefore, comparing the experimental 
values of v 2 with model calculations which do not cor- 
rectly include these fluctuations might not be a priory 
justified. Furthermore, there are effects due to the fi- 
nite width of the centrality bins where also (v 2 } ^ (vz) 71 . 
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FIG. 2: a) Contour plot of the calculated eccentricity in a 
Monte Carlo Glauber model versus the impact parameter, b, 
in Au+Au collisions. The mean value of the eccentricity is 
indicated by the dashed curve, b) The mean eccentricity (e) 
(dashed curve) and the corresponding (e") 1 '" for n = 2, 4 
and 6 (points) versus impact parameter 6. 



The STAR collaboration [HI found a negligible bias on 
the extracted V2 due to these binning effects, except in 
the most central bin. 

To estimate the possible effect of the fluctuations on 
the measured elliptic flow values we calculate the fluctu- 
ations in the initial spatial anisotropy of the created sys- 
tem using a Monte Carlo Glauber model (MCG). This 
anisotropy, which generates the elliptic flow, is given 
by HQ: 



(3) 



where x[ and yl are the coordinates of the constituents 
in the plane perpendicular to the beam and x 1 is in the 
reaction plane (see Fig.^l. Due to the relation between 
the initial spatial anisotropy and the elliptic flow, V2 oc 
e 0j0j^3> fluctuations in e will lead to fluctuations in V2- 
It should be noted that this is only one specific example 
of fluctuations which could contribute to the measured 
value of v% . 

The MCG approach allows for an event-by-event cal- 
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FIG. 3: Eccentricity cumulants calculated for Au+Au colli- 
sions in a Monte Carlo Glauber model compared to (e) as a 
function of centrality (from left to right is from central to pe- 
ripheral collisions, respectively). The symbols are the Monte 
Carlo Glauber results using nucleons, the lines are for con- 
stituent quarks (see text). 



culation of e and therefore the determination of the cor- 
responding higher order moments of the event averaged 
distribution of e. For details on the MCG model and the 
parameters used see Ref. |l9|. 

In Fig. the eccentricity calculated using the MCG 
is plotted versus the impact parameter (b) of an Au+Au 
collision. This figure shows that the fluctuations in e for 
the most peripheral collisions are large. For the most cen- 
tral collisions e can be both positive and negative. This 
leads to an obvious bias when calculating e from (e 2 ) 1 / 2 . 
Figure [2b shows the calculated (e) and the correspond- 
ing (e™) 1 /™. It is clear that there is a bias over the whole 
centrality range, however this bias is the largest for the 
most central and most peripheral collisions. 

Using V2 oc e, we replace in Eq.|21w2 by e, which allows 
us to calculate e{2}, e{4} and e{6}. If the right-hand side 
of Eq. [21 becomes negative, we take the n th root of the ab- 
solute value and multiply this by —1. Figure [3] shows the 
ratio of e{m} and (e), where m = 2,4, 6, versus the colli- 
sion centrality in terms of cross-section. In the standard 
MCG approach for nucleons with a cross-section of 42 mb 
the event by event fluctuations lead to a ratio e{2}/(e) 
which is always larger than 1. This shows that the exper- 
imental determination of V2, using two-particle azimuthal 
correlations in the case of event-by-event fluctuations in 
e, leads to an overestimation of the true elliptic flow. 
This is particularly true for the most central and most 
peripheral events. The inset in Fig.|3shows the same cal- 
culated values on a expanded scale. Even for the central- 
ity region of 20-60% the two-particle correlation method 
overestimates the true elliptic flow by about 10%. In the 
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same figure the results for the higher order cumulant ra- 
tios, e{4}/(e) and e{6}/(e), are shown. The higher order 
cumulants in the case of event by event fluctuations un- 
derestimate the true elliptic flow in the centrality range 
of 0-80%. Above 80% of the cross section also the higher 
order cumulants overestimate the elliptic flow. In the 
mid-central region (20-60%) the higher order cumulants 
underestimate the true elliptic flow by about 10%. If the 
fluctuations are small it can be shown analitically that 
the true value of v-i is given by vi\2} + v 2 {4}/2 |2£|. 

The most central 10% of the cross section show an 
interesting behavior for the different cumulants. The two 
particle correlations always lead to a real value of e{2}, 
however the higher order cumulants can become complex 
because of the combination of two, four and higher multi- 
particle contributions. This is the case for e{4} and e{6} 
in the centrality bin 5-10%. However, in the most central 
bin, 0-5%, e{6} becomes real again. Complex values of 
V2 are usually not reported by experiments. 

To investigate the sensitivity to the magnitude of the 
fluctuations in the eccentricity we followed the work of 
Eremin and Voloshin [2l| and calculate e{m} using va- 
lence quarks as constituents in the MCG. The cross sec- 
tion used for the quarks was set to 6 mb and the num- 
ber of constituent quarks in a gold nucleus to 591. All 
the other parameters were kept the same as in the nu- 
cleon MCG. Due to the larger number of interacting con- 
stituents the fluctuations in central collisions are reduced. 
Indeed, the e fluctuations for the 0-80% of the cross- 
section are reduced as shown by the lines in Fig. 01 

The ratios shown in Fig. |2| cannot be directly com- 
pared to the elliptic flow measurements since {1)2} is not 
an experimental observable. Instead, we calculate the 
ratio e{m}/e{n}. Because of the assumed proportional- 
ity between e and V2, which in principle could depend 
on centrality [TtI ]. this ratio can be compared with the 
measured V2{i r n} / V2{n} . 

Figure^ shows the ratio e{4}/e{2} for the quark and 
nucleon MCG calculations. In both cases this ratio is 
smaller than unity over the whole centrality range, with 
the largest suppression in the case of the nucleon MCG. 
The elliptic flow cumulant measurements from the STAR 
collaboration ^3] are also shown in Fig. The experi- 
mental results are in between the calculated values, and 
are closer to the nucleon (quark) MCG results for periph- 
eral (central) collisions. 

If the observed difference between t>2{2} and ^{4} 
would be due to non-flow effects, it is proposed in 
Rcf. |r| that the squared difference is given by 



■4W = - 

2\ f N 



(4) 



where N is the multiplicity and g%/N is the non-flow con- 
tribution. NA49 9] has indeed observed such a constant 
behavior. 

This squared difference, as measured by the STAR 
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FIG. 4: a) Comparison of the calculated e{4}/e{2} for both 
the nucleon (solid line) and quark (dashed line) MCG to the 
STAR measurements 17]. b) Non-flow or fluctuation contri- 
bution to the azimuthal correlations (see text). 



collaboration is shown in Fig. 0Jd for the cumu- 

lant measurements as a function of centrality. Instead 
of scaling with the multiplicity as done by NA49 we 
scaled the STAR results using the corresponding number 
of wounded nucleons (Nw), which is proportional to the 
multiplicity. These measurements are compared with the 
values obtained from the quark and nucleon MCG (his- 
tograms in Fig.^J)). However, to make this comparison 
the proportionality between V2 and e was assumed to be 
0.16, independent of centrality. It is seen from the figure 
that the centrality dependence of the data is better de- 
scribed by the nucleon MCG than by a constant value as 
expected for non-flow contributions. 

In summary, experimental measurements of elliptic 
flow (^2) might be affected by fluctuations. In this study 
we investigate the effect on the determination of w 2 of 
fluctuations in the spatial anisotropy of the overlap region 
(e) in heavy-ion collisions using a Monte Carlo Glauber 
simulation. For this study we assume a proportional re- 
lationship between e and V2- 

When elliptic flow is measured from two-particle cor- 
relations (v2{2}) we find that fluctuations in e lead to an 
overestimation of its value. On the other hand a cumu- 
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lant analysis of order four (i>2{4}) and six (^2{6}) leads 
to an underestimation of the true value. 

Experimental results on the ratios ^{m} / V2{n} can 
be compared directly to the model and are shown to be 
in approximate agreement. Traditionally, a ratio differ- 
ent from unity is interpreted as being due to non-flow 
effects. However, this study shows that this difference 
can also be explained by event-by-event fluctuations in 
the elliptic flow. In reality it is likely that both non- 
flow and fluctuations affect the measured elliptic flow. It 
remains to be determined which effect dominates. 

Experimentally, this can be addressed by also report- 
ing the ratios V2{m}/v2{n} obtained with the cumulant 
analysis as a function of centrality when they become 
complex because this is a unique feature of fluctuations. 
Event-by-event fluctuations in elliptic flow have been sug- 
gested as an interesting physics observable however 
to address this these more "trivial" causes of fluctuations 
have to be understood first. 
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